A Z2 index of Dirac operator with time reversal symmetry 
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With time reversal symmetry a Dirac operator has vanishing index and Chern number. We show 
that we can nevertheless define a nontrivial Z2 index as well as a corresponding topological invariant 
given by gauge field, which implies that such a Dirac operator is topologically nontrivial. 
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The idea of topological invariants has been successfully 
applied to various fields in physics. In gauge theories 
they are used to classify topological configurations such 
as monopoles and instantons. When couplings with chi- 
ral fermions are introduced, there arise interesting field 
theoretical phenomena such as chiral anomaly and gauge 
anomaly 1] . The chiral anomaly [1] is known to have inti- 
mate relationship with the index theorem [2] which tells 
that the index of the Dirac operator coincide with the 
second Chern number. The gauge anomaly also has topo- 
logical origin, since it is related with the chiral anomaly 
in six dimensions [H- 

In condensed matter physics, it is well-known that the 
plateaus of the quantum Hall effect (QHE) are classi- 
fied by the first Chern number [1, 0] . Recently, a novel 
topological number has been proposed by Kane, Mele, 
and Fu & 0, Q for quantum spin Hall effect (QSHE) 
@i S H Si- is invariant only modulo two, and is of- 
ten called Z2 invariant. Here, time reversal symmetry is 
the key difference between QHE and QSHE. Remarkably, 
QSHE have recently been observed in several experiments 

The formula of the Z2 invariant proposed by Fu and 
Kane [8] is, roughly speaking, "half" the first Chern num- 
ber. Therefore, it is very useful IJ] for numerical calcula- 
tions if we utilize the techniques of cornputing the Chern 
number in the lattice gauge theories [l5|. Besides such 
practical applications, it is of fundamental importance, 
since it could be a topological invariant in "mod 2 index 
theorem" fl6l|. Therefore, if we find a corresponding an- 
alytical invariant, we can obtain a simple formula of a Z2 
index theorem for a (pseudo-real) Dirac operator. 

In this paper, we investigate analytical and topological 
invariants associated with a Dirac operator with time 
reversal invariance. We first study its spectral properties 
in Euclidean space and define a Z2 index of the Dirac 
operator. We then propose a topological invariant which 
is a generalization of the Fu-Kane formula, and infer that 
it coincides with the Z2 index. 

We begin by recalling the time reversal transforma- 
tion of Dirac fermions in d = 2?i + 2 {n = 0, 1, • • • ) 
dimensional Minkowski space-time. It is defined by 
ip{t,x) — > Til![~t,x), where T is an anti-unitary op- 
erator, T = r-yFc/C, with F-y being a product of some 
7 matrices, Fq a generator of a gauge group G, and 



/C the operator of taking complex conjugate. For the 
Lagrangian density C{t,x) — ■ip{t,x)ip{t,x)tlj{t,x) to 
transform as C{t,x) — > C{—t,x) under time reversal, we 
seeTip{t,x)T~^ — ip{—t,x), from which it follows that 



TA^{t,x)T-^ ^A''{-t,x) 



(1) 



where the metric is g^^ = diag(l, —1, • • ■ , —1). The F5 
matrix anti-commuting with the Dirac operator is given 
by F5 i<i/2-i^o^i . . .^d-i^ rpj^jg definition directly 

leads to TEgT^i = -F5 for d = 4n + 2 and = +F5 
for d = 4n -I- 4. 

T has the following two possibilities: — ±1, de- 
pending on F^ = ±1 and Fq = ±1. The former is deter- 
mined solely by the space-time dimension d: In d = 2, 
for example, we can choose 7*^ = and 7^ 
these are imaginary, IC{j'^ ,j^)IC^^ — (— 7^^ 



that F, 
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we see 

-17^ = a^, and therefore, F^ = 1. In gen- 
eral, we have F^ = 1 for d = 0, 2 -I- 8n, and F^ = — 1 for 
d = 4,6 + 8n. 

To discuss the index of the Dirac operator with time 
reversal invariance, we switch from Minkowski space to 
Euclidean space. It should be noted that a Euclidean 
version of the time reversal transformation is not so ob- 
vious, since it includes the operator /C. Here, we define 
Euclidean space by rotating all the spatial coordinates 
(j = 1, • • • , c? — 1) onto the imaginary axes via = iy^ , 
whereas x^ — y'^. As we shall see, this enables us to relate 
a Z2 index of the Dirac operator with a topological invari- 
ant [17]. The metric becomes in this case g^^ = 5^i,. Cor- 



respondingly, the 7 matrices are introduced via 7-' 



and 7 



which become hermitian 7 



27-' 



and 

the gauge potential A^{y) via Aj{x) — ~iAj{y) and 
Ai^{x) — Ad{y)- Then, the Dirac operator can be de- 



noted as ip{x) — ip{y) = ij'^idy^ ~ 
regard as hermitian {ip{y)y — ip{y). 
tion law under time reversal becomes 



iA{y)) which we 
The transforma- 



TA^{y)T-^=A^{-y), 



(2) 



which follows from the same transformation law ([T]) but 
with the Euclidean metric mentioned above. Therefore, 



2 



the Dirac operator transforms as 



(3) 



Note that this transformation is for flat space: If we con- 
sider curved space, it should be modified suitably, as we 
shall see. Other key symmetry is chiral symmetry de- 
scribed by 



(4) 



The Z2 index discussed in this paper is involved in the 
case with both of the conditions 



= -1, 



(5a) 
(5b) 



fulfilled. The former ensures that the eigenstates are al- 
ways doubly-degenerate, which is referred to as Kramers 
doublet. The latter claims that the zero- mode Kramers 
doublet have opposite chiralities. These conditions give 
some constraints: Eq. (|5bp is valid only ind = An+2, and 
Eq. dSaj) imposes = -1 for = 8n + 2 and T^ = +1 
for d 8n -f- 6. Typical example in the former case is 
Eg — In ® = J2, whereas in the latter case, a con- 
venient but nontrivial choice may be Eg = Iat = Ji 
The transformation law of gauge potentials is thus 
defined by 



J2A;{y)J^ 



for d = 8n 



(6) 



Eor a time reversal invariant Dirac operator discussed 
so far, we shall define a Z2 index. Let ipkiu) be an eigen- 
states of ip: 

ip{y)(pk{y) = ek(pk{y)- 

Then, Eq. ([3]) ensures that (pKk{y) = Tvk{—y) is also an 
eigenstate of ip{y) with the same eigenvalue Sk- Here, 
the condition ([5a|) plays a vital role in the orthogonal- 
ity between ipk and (pKk- It thus turns out that all 
eigenstates are doubly-degenerate, called Kramers dou- 
blets as mentioned above, which we denote as ^k{y) = 
{ipk{y), ^Kk{y))- The spectrum is illustrated in Fig. [TJ 

Let us concentrate on the zero-mode eigenstates, $o,a 
(a = I,-- - ,N-k). Since the Dirac operator anti- 
commutes with E5, the zero- modes can be chosen to be 
eigenstates of the chirality. Suppose T5ipo.a — +1^30.0. 
Then, we see T^ipKo.a = —'fiKo.a because of Eq. (j5bp . 
Namely, at the zero energy, each Kramers doublet is com- 
posed of two states with opposite chiralities. Even when 
there are some doublets at the zero-energy, the number 
of states with positive chirality is the same as the number 
of states with negative chirality: 

ind?.^ = dimkeri^_|_ — dimkeri^_ — 0, 



r 



FIG. 1; Schematic illustration of the spectrum of the Dirac 
operator in the case TVk ~ 2 (two zero-mode doublets) on 
a compact manifold M. The nonzero-mode quartet can be 
obtained by the operation of T and /or Fs. The chirality ± is 
shown for the zero-mode eigenstates. 



where ip± = ipP± with P± = (1 ± E5)/2. The index 
of the present Dirac operator is thus trivial. Neverthe- 
less, the time reversal invariance ([3]), if combined with 
chiral symmetry ([4]), gives an interesting invariant. Chi- 
ral symmetry (jl]) tells that if ^k is an eigen-doublet with 
the energy e^, the state defined by ^-k = ^h'^k is also 
an eigen-doublet with the opposite energy — e^. There- 
fore, nonzero-mode states form a quartet in this sense. 
Suppose that we have just one Kramers doublet at the 
zero energy. Then, it turns out that this doublet is stable 
against perturbations with time reversal and chiral sym- 
metries, since these two states cannot move to nonzero 
energies without two more states in order to ensure both 
the symmetries. On the other hand, if there are two dou- 
bles at the zero-energy, they are not obliged to stay there: 
Small perturbations enable two of them to move to the 
positive energies and the other two to move to the op- 
posite negative energies. In more general, we can claim 
that evenness or oddness of the number of the zero-mode 
Kramers doublets is an analytic invariant, from which we 
define a Z2 index of the Dirac operator with time reversal 
symmetry. 



ind-|_ ip — dimkeri_^+ mod 2. 



(7) 



Next, we define a topological invariant given by the 
gauge field. To this end, we must first specify the mani- 
fold compatible with the condition Q . Let M be a com- 
pact manifold without boundary. We assume that it can 
be divided into two M± such that \i y € M_|_, — y G M_ 
except for the time reversal invariant points yj — ^yj 
{j = 1, • • • , Mnv)- The number TVjnv of such points de- 
pends on M. For example, the two-sphere S^ has two 
invariant points, whereas the two-torus T^ has four in- 
variant points, as illustrated in Fig. [2l 

To define a topological invariant, it is convenient to 
define a gauge potential one-form A — —iA^j^Ay^ and cor- 
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FIG. 2; An example of two-dimensional manifold M. de- 
notes the time reversal invariant points. The square rep- 
resents if the boundary is regarded as one point. In 
this case, time reversal invariant points are just two, j/i and 
2/2 (= 2/3=2/4)- K the two parallel boundaries are pasted and 
the periodic boundary conditions is imposed on each direc- 
tion, the same square now denotes T^, which has four time 
reversal invariant points. 



responding field strength two- form T = + . Time 
reversal invariance ([2]) or ([6]) tells that the Chern number 
Cd (d = 4n + 2 with n = 0, 1, • ■ • ) is vanishing: 



M 



where the numerical factor is ,y\^n = i"^ /\n\{2n)"'\. This is 
consistent with the spectral property of the Dirac opera- 
tor whose index is zero, as discussed above. This feature 
is quite similar to the QSHE: The first Chern number 
associated with the Berry phase in the Brillouin zone 
(T^) vanishes due to time reversal symmetry. Neverthe- 
less, the QSHE phase is the topologically nontrivial phase 
which can be described by the Z2 number, as shown by 
Kane et. al. 0, Q ■ Motivated by their work, we propose 
that the Z2 index ([7]) in d = 4n -|- 2 is equivalent to 



Dd^^d 



d/2 



M4 



am. 



(8) 



mod 2. Here, LOd-i is the Chern-Simons {d — l)-form 
which obeys ir T'^/'^ = AuJd-i{A) It should be noted 
that the formula ([5]) have meaning under the condition 
This is the reason_we have adopted an unconven- 
tional Wick rotation 
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We must first examine the gauge-dependence of Dd- 
Let Ag = g^^Ag + g~^dg be the gauge-transform of 
A. If the time reversal invariance ([6]) is enforced on Ag, 
g should obey g{—y) = J2g*(y) for ci = 8n -I- 2 or 
g{-y) = Jig*{y)Ji^ for d^8n + 6. We will refer to this 
condition as the time reversal constraint on the gauge 
transformation. Let Ad_i[(/] = £)d[.4g] — -Dd[-4] be the 
gauge-dependence of Dd- Note 

ujd-i{Ag)~uJd-i{A) = ojd-iig^^dg) + dad-2, 



where ad-2 is a (d — 2)-form ^lo'l, which leads to 



(-1) 



n+l- 



(2n) 



(27r)2"+i (4n+ 1)! 



tr(g-Mg) 



4n+l 



where d = An + 2. Let us estimate the above in the 
case M = S'' (9M+ = S'^-^) for simphcity. Notice that 
generic U(2A^) gauge transformation g can be decom- 



posed into U(l)xSU(2iV) such that g{y) 



oi<t>(y)n 



g{y) 



where det g ~ 1. The time reversal constraint tells that 
(j){-y) = -0(2/) mod 27r. In d = 2 (i.e., n = 0), Ai[5] 
is given by this U(l) part, Ai[(/] = —N/ir § where 
line integral is over S^, namely, the equator of S^. This 
gives manifestly an even integer. On the other hand, 
in higher dimensions, contribution from U(l) vanishes 
and only the non-Abelian sector g enters into Ad_i; 
Ac;-i[(/] = /S.d~i[g\. The time reversal constraint tells 
that at the time reversal invariant points yj (j' = 1, 2 on 
S''"^), giVj) = bh/j) e Sp(A^) for d = 8n + 2, whereas 
h{yj) e 0(2iV) [20] for d ^ 8n + 6. First, let us consider 
the former case. Assume that g takes ^0 ^ Sp(2A^) at a 
certain y, g(y) = go. Then, y cannot be yj, and the time 
reversal constraint ensures that at g{—y) — go- There are 
thus even number of points on S''"^ which are mapped to 
go- It turns out that the degree of the map g is even, im- 
plying that the winding number, Ad_i [5] , is even. On the 
other hand, if one cannot find any ^0 ^ Sp(iV) on S'^"^, 
namely, if g{y) G Sp(iV) for all y, Ad-i[g] = 0. We thus 
conclude that Ad-i[g] is an even integer for d = 8n -I- 2. 
The case d = 8n -I- 6 (n = 0, 1, • • • ) is likewise. 

If the gauge potential can be smooth on the whole 
M+, Dd should be zero, which is a trivial element of Z2. 
Now we shall show that there exist not only such a trivial 
element but also a nontrivial element indeed. We assume 
that the u(2A^) gauge potential is 2 x 2 block-diagonal. 
Then, the time reversal invariance ([6]) requires that the 
upper and lower u(A'^) sector of the gauge potential is not 
independent, given generically by the form 



AM 



(9) 



where 0^(2/) denotes a u(A^) gauge potential. In this 
block-diagonal case, since the upper and lower sectors 
are decoupled, the Z2 index and Dd can be separately 
computed such that md+ip = ind+ip^ + ind+ipi and 
Dd — Dd] + Ddi , where arrows mean that only the up- 
per (1) or lower (|) gauge potential is taken into account. 
Assume that the gauge potential in the upper sector is 
nontrivial, which yields a nonzero index and Chern num- 
ber. This means that a^{y) cannot be smooth over S'^: 
For simplicity, assume that we have two kinds of gauge 
such that a'jt^\y) is regular in M±, collectin g; al l singu- 



larities in Mzp. The ordinary index theorem [2l| claims 
that ind ip] = — N_, where N± is the number of 
the zero- mode with the chirality ±. Now let us choose 
ajr^ (y) as the upper gauge potential. Then, Eq. ([5]) gives 
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D 



ind 



N^, since in this case, \y) is 



regular in M_ and therefore, 

JdM+ JdM- JM- 

holds if the lower gauge potential is neglected. Next, let 
us switch to the case with lower gauge potential only, 
which should be a^~'*(— y). Since this gauge potential is 
regular in M+, it never contribute to Eq. ([5]); D^j, = 0. 
It thus turns out that = — N_ holds for the full 
gauge potential Af^ in Eq. © . On the other hand, since 
ipl has jVzp zero-mode with chirality ±, i.e., indilpi = 
iV_ - 7V+ dH, we reach ind+ ip ^ N+ + N^. Therefore, 
we conclude that ind+ ip — Dd mod 2. 

Let us now take into account off-diagonal elements of 
the gauge potential. The singularities in the upper gauge 
potential may move but stay in M+ if off-diagonal ele- 
ments are small enough. Even if one of them moves into 
M_ , its partner in the lower gauge potential in M_ moves 
into M+. This is due to Eq. ([6|) which ensures that if A^j. 
has a singularity at j/, an opposite singularity appear at 
—y. Therefore, Dd can change only by 2. On the other 
hand, along the change of the gauge potential, the spec- 
trum of the Dirac operator flows, and a nonzero-mode 
quartet can be two zero-mode doublets and vice versa, 
which result in the change of Z2 index also by 2. From 
the point of view of such moving singularities, the mod 2 
gauge-dependence of Dd can be understood likewise. It 
thus turns out that Dd and the Z2 index change by two 
and therefore coincide mod 2. 

Finally, we shall exemplify a Dirac operator with non- 
trivial Z2 index in d = 2. Let us consider a Dirac operator 
on with magnetic monopole background fields 22l |. 



ip{e,^ 



ia^ I da 



— cot 6 — iAf 
2 



(90 - iAf/,) 



where < ^ < tt and — tt < (j) < tt are polar coordinates, 
Ae and A^ are u(2) gauge potentials of the type ([9]), and 
the cotangent term is due to the spin connection [j^ . 
This Dirac operator has time reversal symmetry 

Tip{e, (l))T-^ ^ ip{TT - e, -0), 

where T = a^ir^lC. Two time reversal invariant points 
are (0,0) = (7r/2,0) and (7r/2, tt). For the upper gauge 
potential we have two well-known possibilities. 



.(±) 



(0) 



-(±1 



COSf 



and ag — 0, where a^^^ is the charge-m monopole poten- 
tial with a singularity at the south and the north pole, 
respectively. From Eq. ^ it follows that the lower po- 
tential should be a^'^'*(7r — 0) ^ ~a'^\o), telling that it 
denotes a monopole with the opposite charge and with 
the singularity at the opposite pole. 



Assume m > and choose a^~^ as the upper gauge 
potential. Then, ip-[ gives just m zero-modes with chi- 
rality whereas ipi gives the same m zero-modes but 
with chirality ^. On the other hand, we see D^t = ^ 
and Ddi = 0. Therefore, for this decoupled model and 
the present gauge-fixing, the Z2 index and D2 coincide, 
ind+ = m = Z?2- However, as discussed, these can 
change by 2 by gauge transformations and/or deforma- 
tion of the gauge potential, and generically coincide mod- 
ulo 2. 
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